Thermodynamics of a dry atmosphere at different 
surface exchange rates and rotation speeds 



Salvatore Pascale a , Francesco Ragone a , Valerio Lucarini a ' b , Yixiong Wang c 

a KlimaCampus, Meteorologisches Institut, Universitdt Hamburg, Hamburg, Germany 
b Department of Mathematics and Statistics, University of Reading, Reading, U 
c University of Oxford, Clarendon Laboratory, Parks Road, Oxford, UK 



Abstract 

We study the combined effect of the rotation speed Q and of the surface ex- 
change rate - quantified by a surface turbulent relaxation timescale r - on the 
dissipative properties of an Earth-like dry atmosphere. The rotation speed Q 
is varied between one tenth and eight times that of the Earth Qe ~ 7.29- 10 -5 
rad -1 and r from 45 minutes to 500 days. We study the circulation regimes 
induced by such parametric variations through two key dimensionless pa- 
rameters, the thermal Rossby number TZo and the frictional dimensionless 
number J 7 /. An extensive analysis is performed by using nonequilibrium 
thermodynamics diagnostic tools such as material entropy production, effi- 
ciency, meridional heat transport and kinetic energy dissipation. The thermal 
dissipation associated with the sensible heat flux is found to depend mainly 
on the surface properties and to be almost independent from the rotation 
rate, whereas the dissipation of kinetic energy depends in a nontrivial way 
on both. Slowly rotating, axisymmetric circulations (TZo > 1) have the high- 
est mechanical dissipation when the surface drag is strong (J 7 / ps 10~ 3 ), but 
the highest efficiency for T$ ~ 10. For 0.01 < TZo < 1 the peak is reached for 
Tf ~ 10 3 (r ~ 3 d), corresponding to the maximum activity of the baroclinic 
eddies, the maximum meridional heat transport and the highest efficiency. 
At high rotation rates (TZo < 10~ 2 ) there is a dramatic drop in the intensity of 
the atmospheric energy cycle and in the meridional heat transport as the at- 
mosphere tends towards the radiative-convective equilibrium profile. When r 
is interpreted as an internal parameter, our results also confirm the vagueness 
of the Maximum Entropy Production Principle, since its applicability seems 
to be dependent on both the dissipative functions and the dynamical regime. 
This study suggests the effectiveness of using fundamental nonequilibrium 
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thermodynamics to investigate the properties of planetary atmospheres. 
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In the last two decades, more than 700 planets outside the solar sys- 



tem (exoplanets) have been discovered (Dvorak, 2008; Bonfils and coauthors 



2012), and the Kepler Space Telescope has located over 2,000 exoplanet can- 
didates. The study of exoplanets and their climates is in its early stage and 
it is quickly developing. Observational data are still poor and difficult to 



obtain, particularly for those planets - super-Earths (Valencia et al. 2007) 



- that might be capable of sustaining liquid water and thus potentially suit- 
able for life. Nevertheless, the discovery of exoplanets is pushing the limits 
of planetary sciences, so far concerned with the study of the Earth and solar 



system planets, towards very diverse climate scenarios (Heng, 2012 A. Bur- 



rows 


1997; 


Heng et al. 


2011a; 



2009). Exoplanet climates 
are in general characterized by a wider range of physical - e.g. atmospheric 
composition, rotation rate, dimension, surface - and orbital - e.g. obliquity, 
eccentricity, distance from the parental star, spectral type of the parental 
star, phase locking - parameters than planets belonging to the solar system 
(Williams and Pollard 2002), and therefore span an even greater diversity of 
circulation regimes. Following Read (2011) and Wang (2012), the variety of 
circulation regimes may be reduced by adopting the fluid-dynamical classical 
method of similarity, i.e. by defining a set of dimensionless numbers that 
fully characterise the climate states. Two climate states that share the same 
set of dimensionless numbers are then dynamically equivalent, so that the 
statistical properties of one can be mapped onto those of the other. Obvi- 
ously the number of parameters is fairly large, and one of the main objectives 
of planetary science is to understand what is the minimal number of parame- 



ters needed in order to characterise effectively the general circulations (Read 



2011 Showman et al. 2010). 



In this study we are going to focus on the rotation rate and on the sur- 
face turbulent exchange rate of an Earth-like dry climate. Over the last three 
decades, the effect of the planetary rotation has been investigated in some de- 



tail with the aid of general circulation models ( 


Hunl 


1979 


Williams 


1988alb 


Navarra and Boccaletti 


2002; 


Genio and Suozzo 


1987 


; 


Gleisler et al. 


1983 


Read 


2011 


Vallis and Farneti 


2009). Such studies revealed a clear trend in 
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the way the general circulation of an Earth-like planet depends on changes in 
Q. The terrestrial circulation features an axially symmetric Hadley regime 
at low latitudes and quasi-geostrophic regime at higher latitudes charac- 
terized by baroclinic eddies, with both contributing to the meridional heat 
transport. The most direct effect of O is on the size of the baroclinic distur- 
bances, which decreases like the Rossby deformation radius as I/O, and on 
the latitudinal extent of the Hadley cell which also decreases with the rota- 
tion rate according to I/O, as shown by |Held and Hou (1980). Simulations 
of slowly rotating Earth-like planets (typically at rotation rates smaller than 
half of the Earth's) show the presence of one axially symmetric direct cell 
(Hadley cell) in each hemisphere, which becomes the dominant mechanism 
for meridional heat transport, whereas the quasi-geostrophic baroclinic cir- 
culation ceases to exist in the middle and high latitudes. On the contrary, at 
high rotation rates we observe the emergence of multi-cell structures in the 
meridional circulation (Williams, 1988a). 

The dynamical effects of the solid lower boundary of terrestrial planets are 
among the specific research issues which have not been fully addressed yet, as 
discussed in Showman et al. (2010). The characteristics of the surface have 



been recognised as a key factor in shaping the atmospheric circulation ( James 



1994 James and Gray, 1986 Read, 2011), although this topic has received 



less attention than the rotation rate. The surface of a terrestrial planet has 
a considerable effect on the turbulent flow within the planetary boundary 
layer because of its roughness (Arya, 1988). This affects the exchange of 



momentum and energy between the surface and the free atmosphere and 



consequently the whole general circulation. It has been shown (James and 



Gray 1986 James 1987 Kleidon et al. 2003[ ) that the reduction of the 
surface drag leads to strong horizontal barotropic (i.e. independent of height) 
shears in the zonal mean flow. By using a two-level quasi-geostrophic model, 



James (1987) showed that the growth rate of the most unstable baroclinic 



modes is reduced considerably by the strong horizontal wind shears. This is 
related to the general fact that the linearised baroclinic instability equations 
obey the Squire's theorem (Kundu and Cohen, 2004). 

However, to the authors' knowledge, the effect of surface drag has not 
been systematically investigated so far at rotation rates which are different 
from the Earth's. In this study we investigate the combined effect of rotation 
speed and surface roughness on the thermodynamics of the system. Ther- 
modynamics provides a way to characterize in a concise way a climatic state, 
bringing together comprehensive but minimal physical information. This is- 
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sue is closely related to the definition of robust observables (Lucarini, 2012). 



The atmosphere of a planet is a nonequilibrium system (Gallavotti, 2006 



DeGroot and Mazur, 1984 Kleidon, 2009), and its general circulation redis- 



tributes energy in order to compensate for the radiative differential heating 
between hot and cold regions - typically low- and high-latitude regions in 
the case of e.g. Earth, Mars, Titan (the other way around for planets with 
large obliquity such as Uranus). The atmospheric circulation therefore is 
fuelled by the conversion of available potential energy due to large tempera- 
ture gradients into kinetic energy. The atmosphere, in other terms, produces 



mechanical work, acting as a heat engine (Lorenz, 1967; Peixoto et al. 1991 



Johnson , 


2000 


Lucarini , 


2009) 



In the Earth system the disequilibrium drives a variety of irreversible pro- 
cesses, from air motion with the subsequent frictional dissipation to chemical 
reactions. The irreversibility of climatic processes is quantified by the ma- 



terial entropy production (Goody, 2000 Kleidon and Lorenz, 2005; Kleidon 



2009). The interest in studying climate material entropy production largely 



stemmed from the proposal of the maximum entropy production principle 
(MEPP) by Paltridge QPaltridge[ |1975[ |1978| |2001[ ), which suggests that the 
climate adjusts in such a way as to maximize the material entropy produc- 
tion. Whereas its theoretical foundations are still unclear (Dewar, 2005 



Grinstein and Linsker, 2007 Goody, 2007), such a conjecture has also been 



proposed as a way to estimate the meridional heat transport of other planets 



too, such as Mars and Titan (Lorenz et al. 2001 Jupp and Cox, 2010) and 



has stimulated the re-examination of climatic dissipative processes (Peixoto 



eTaLl |i"99T| |Goody[ [20001 IPauluis and Heldj [2002iHb| |Kleidon and Loreni 



2005; Fraedrich and Lunkeit 



2008 



Pascale et al. , 2011a). 



Recently, a link has been found among the Carnot efficiency, the intensity 
of the Lorenz energy cycle, the material entropy production, and the degree of 
irreversibility of the climate system ( Lucarini , 2009 ) - namely, the efficiency 
of the equivalent thermal machine also sets the proportionality between the 
internal entropy fluctuation of the system and the lower bound to entropy 
production by the fluid compatible with the second law of thermodynamics. 
Such a bound is basically given by the entropy produced by the dissipation 
of the mechanical energy, whereas the excess of entropy production is due 
to the turbulent transport of heat down the gradient of the temperature 



field. Such an approach has been insightfully applied for example in Lucarini 



et al. (2010) and Boschi et al. (2012) for studying the bistability of the Earth 



system associated with the sea-ice feedback. 
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In this study we extend such an analysis based on the use of fundamental 
nonequilibrium thermodynamics to investigate the properties of planetary 
atmospheres and consider two other fundamental parameters in shaping the 
general circulation of the atmosphere: the planetary rotation speed and sur- 
face exchange rate. 

The structure of the paper is as follows. In Section [T] we will shortly 
discuss the dimensionless parameters approach and put our simulations in 
such a context. In Section [2] the model and the experimental setup is pre- 
sented, and the basic thermodynamic tools recalled in Section [3j Results are 
presented in Section [4], and finally, in Section |5j the main conclusions are 
summarized. 



1. Parametric range of general circulations and dimensionless num- 
bers 

The role of the rotation rate in planetary circulations has been inves- 
tigated first in laboratory experiments with the thermally driven rotating 



2001 Wordsworth et al. 2008 Hide, 2010), which is capable of capturing 



1975 


Read et al. 


1998; 


Read 



some of the basic features (horizontal differential heating and the planetary 
rotation) of the atmosphere's general circulation. The system consists of a 
fluid confined between coaxial cylinders maintained at two different temper- 
atures and rotating at an angular velocity Q. When the basic parameters Q 
and AT (temperature difference between the inner and outer cylinder) are 
varied in a controlled way, a wide variety of flow patterns are observed. Such 
patterns show some clear trends and transitions between different regimes if 
results are grouped with respect to two dimensionless parameters, the ther- 
mal Rossby number: 

gaDAT 



and the Taylor number: 



Tlo 



Ta 



2 r 5 



4Q 2 L 
v 2 D 



(1) 



(2) 



in which L is the channel width, D its depth, v the kinematic viscosity of the 
fluid, a its volumetric expansion coefficient, and g the gravitational acceler- 
ation. The parameter Ta measures the strength of the viscous dissipation 



with respect to the Coriolis force and turns out to be (Read, 2011): 



Ta oc 41TT, 



2_2 



(3) 
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since the viscous timescale is r„ = L 2 jv and L/D ~ 1. On the other hand, 
IZo is regarded as a measure of the strength of buoyancy relative to Coriolis 
forces. 



Read (2011) has proposed that similar dimensionless numbers can be 
defined for atmospheric circulations too, so that the diversity of planetary 
circulations observed in our solar system and eventually on exoplanets may 
be more effectively understood and classified. The analogous thermal Rossby 



number (eq. (fij)) is defined as (Read, 2011): 



RA9 h 

where a is the planet's radius, R the specific gas constant and A9h the hor- 
izontal (potential) temperature contrast between equator and poles. A dif- 
ference between the definitions in eq. ([I]) and eq. Q is that A9h is not fixed 
externally but rather determined by the circulation itself. In the following we 
will take A9h = A9hE, where 9hE is the radiative-convective equilibrium po- 
tential temperature, as this is externally determined by the incoming stellar 
radiative energy and thus a more objective quantity to describe the horizontal 
differential driver for the circulation. A dimensionless parameter analogous 
to the thermal Rossby number can be defined starting from ^ as: 

Tj = (5) 

in which 77 is the typical timescale for kinetic energy dissipation. We will 
focus on this number in the present study. For planets with a solid core, r/ is 
the surface drag timescale and is in general determined by the characteristics 
of the surface. 



2. Model and experimental setup 

2.1. The Planet Simulator 



Experiments are performed with the Planet Simulator (PlaSim, Fraedrich 



et al. (2005)). This is motivated by the fact that our study requires us to ex- 



tensively alter parameters of the climate system, as well as producing many 
simulations. Therefore we need a model that is flexible and fast to run rather 
then a state-of-the-art model encompassing as many Earth system processes 
as possible. PlaSim is a general circulation model of intermediate complexity. 
The three-dimensional equations are solved using a spectral transfer method 
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(Eliasen et al. 1970 Orszag, 1970). Parameterisations for unresolved pro- 



cesses deal with longwave (Sasamori, 1968) and shortwave (Lacis and Hansen 



1974) radiation, interactive clouds (Stephens, 1978 Stephens et al. , 1982 



Slingo and Slingo, 1991), moist (Kuo, 1965, 1974) and dry convection, large 



scale precipitation, boundary layer fluxes and vertical and horizontal dif- 
fusion (Louis 1979; Louis et al. , 1981 Laursen and Eliasen, 1989). In the 



basic configuration, the ocean is represented as a mixed-layer (swamp) ocean, 
which includes a zero-dimensional thermodynamic sea ice model. Recently, 
the possibility of running the model in a fully coupled atmosphere-ocean 
setting has been implemented (Dahms et al. 2012). The model is freely 



available at [http: / /www.mi.uni-hamburg.de/plasim| Moreover, a tested en- 

thus making it 



2008) 



tropy diagnostic is available (Fraedrich and Lunkeit 
well suited for this work. 

All experiments are performed at T21 spectral resolution (approximately 
5.6° x 5.6°) and with five vertical levels. The resolution T21 inevitably leads 
to an underestimation of the strength of the Lorenz energy cycle and of the 
material entropy production (e.g. of a factor 0.7 in Kleidon et al. 2003), as 



the smallest baroclinic eddies are not fully resolved, and it has been chosen 
in order to have a fast running model. In order to simplify any problem we 
choose an idealized topographically flat setup and consider a dry atmosphere. 
In a subsequent study we shall also consider the impact of the hydrological 
cycle and its feedbacks. In order to assess the effect of the climate's heat 
capacity, we perform two sets of experiments where we consider a high and a 
low value for the heat capacity of the surface. In the first case (HIGHHC), the 
model has a heat capacity C ~ 2 x 10 8 Jm~ 2 K _1 , which is a typical value 
for a 50-meter deep ocean mixed layer. In the second set of experiments 
(LOWHC) the heat capacity is C « 2.5 x 10 7 Jm~ 2 K -1 , which is typical 



for a rocky planet (Wordsworth et al. 2011). For both set of runs we take a 



surface albedo equal to 0.2, also typical for terrestrial planets (Wordsworth 



et al., 2011). 



2.2. The strength of the turbulent surface exchanges 

In order to have a wide and controlled variation in Tf (Eq. [5]), we simplify 
the representation of the surface fluxes. In PlaSim the temperature tendency 
of the first atmospheric layer (of thickness dz) due to the turbulent sensible 
heat flux Ft = Ch\u\(T s — jT) is computed as: 



dT 
~dt 



pcpdz 



C h \u\ 

dz 



{T s -lT) 



Ts-jT 

Th(*,t) ' 



(6) 
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where Th has time dimension and is computed at each grid-point x at time 
t following the Louis' scheme (Louis 1979 Lunkeit et al. 2010). The value 



of Th in a standard run is a function of space and time but remains of the 
same order of magnitude. Since we are interested in variations of orders of 
magnitude in r h , we substitute the locally computed 77, with a fixed (in space 
and time) time scale as: 



dT _ jT- T s 
dt T h 

Similarly, for the surface stress we have: 

du 

~dt ' 



u 



(7) 



(8) 



Generally the drag and heat transfer coefficients Cd and Ch - and therefore 
the time constants T m and Th - have similar magnitude. This is particularly 
true in the case of neutral flows, for which Co = Ch is indeed a very good 
approximation (Arya, 1988; Louis, 1979). For non-neutral flows, Ch and Cd 



are different but still of the same order of magnitude, as can be seen in Fig. 
11.6 of Arya (1988). On the base of this and since in this study we are going 



to explore a wide parametric range, we assume for the sake of simplicity: 



T„ 



Th 



T. 



(9) 



Experiments are performed for Vl* = Vl/VL E = 1/10, 1/5, 1/2, 1, 2, 4, 8, 
where Qe = 7.29 • 10 -5 rad -1 is the Earth rotation rate. For each value 
of tt* we run the model with r = 2700, 3600, 10800, 21600, 43200, 86400, 
259200, 864000, 2592000, 8640000 s, that is from 45 min (model timestep for 
Q/Qe < 1 ) to 500 days. Simulations with very large r are representative of 



an atmosphere with no solid lower boundary (James, 1994) and in practice 



used in simulations when general circulation models developed for the Earth 
are adapted to simulate Hot Jupiters (e.g. Menou and Rauscher 2009 Heng 



et al., 2011b) 



Let us note that as Q increases, the typical size of the baroclinic distur- 



bances L c decreases as (Eady, 1949) 



2AirL Rl 



(10) 



with the Rossby deformation radius Lr = NH/f (James, 1994 Williams 



1988a) (N the buoyancy frequency, H the height scale and / = 20, simp 
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the Coriolis parameter). In particular, we have to bear in mind that L c ~ 
6000 km for Q* = 1 and L n ~ 1500 km for fi* = 4 QWiUiams[ |1988a| ), thus 



indicating that L c < 750 km for fi* > 8. At T21 resolution (~ 500 km) 
therefore the model does not resolve anymore baroclinic eddies for Q* > 8. 
Because of this we have not performed simulations at Q* > 8. 

It is useful now to quantify the experiments within the parametric space 
(TZo^f). In the estimates of IZo we have taken A&h ~ 60 K which cor- 
responds approximately to the potential temperature difference near the 
surface of the radiative solution for such a dry atmosphere. The range of 
variation of the dimensionless parameters TZo and JFf is shown in Tab. [TJ 
We see that TZo varies approximately over four orders of magnitude, whereas 
J-f varies over thirteen orders of magnitude, thus extending considerably the 
parametric range for T$ considered by Williams (1988a. b) and Read (2011) 
(4 orders of magnitude for both parameters). Typical values for Earth, Mars 
and Titan of IZo and T s a re 6 x 10" 2 , 2 x 10 _1 , 18 and 1.6 x 10 4 , 44, 7.5 x 10 4 
respectively (Read, 2011). 



3. Thermodynamic diagnostics 

The general circulation is the result of the conversion of the available 
potential energy generated by radiative differential heating into mechanical 



work (winds), as first shown by Lorenz (1955, 1960, 1967). For an atmosphere 
in a statistical steady state, the rate of generation of available potential en- 
ergy, G, the rate of conversion into kinetic energy, W, and the rate of dissi- 
pation of kinetic energy through the turbulent cascade (and ultimately via 
viscous dissipation), D, have to be equal when averaged over long time peri- 
ods (e.g. a year or longer), G = W = D ((•) denotes the time mean). They 
are therefore equivalent ways of measuring the strength of the atmosphere's 



general circulation (Peixoto and Oort, 1992). 



The energy cycle introduced by Lorenz has been set onto a thermody- 
namic framework through the consideration of the effective Carnot engine 



describing the ability of the atmosphere to perform work (Johnson, 2000 



Lucarini, 2009). The atmosphere is seen as a heat engine which generates 



mechanical work at average rate W from the differential heating due to ra- 
diative and material (e.g. latent heat release) diabatic processes. If Q + and 
Q~ are the local positive and negative diabatic heating rate (i.e. Q + = Q 
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where Q > and Q + = where Q < and similarly for Q ) with 



± 



Q ± pdV, 



(11) 



we have that $~ 
as: 



+ <3>~ = W > 0. Moreover, one can define an efficiency 77 



(12) 



which gives us an indication about the capability of the general circulation 
of generating kinetic energy given the net heating input $ + . From Eq. (12) 
it follows that 

W = rjW (13) 



in full analogy with the definition of efficiency of a heat engine (Fermi 1956[ ). 
Such a quantity has been proved to be particularly relevant in marking the 
climatic shifts between the present day climates and the Snowball Earth 



(Lucarini et al. 2010 Boschi et al., 2012) 



Dissipation, and therefore irreversibility, is ubiquitous in planetary at- 
mospheres and, more generally, in nonequilibrium steady state systems. The 
kinetic energy of the atmospheric flow is ultimately transferred through a tur- 
bulent cascade to smaller scales where it is then dissipated into heat by fric- 
tion due to viscosity. Thermal dissipation due to sensible heat fluxes between 
the surface and lower atmosphere is another irreversible process which may 
take place in planetary atmospheres. Planets whose atmospheres allow phase 
transitions of one or more of their chemical substances (e.g. water on Earth 
or methane on Titan) also experience further irreversible processes as evap- 



oration/condensation and diffusion (Goody, 2000 Pauluis and Held, 2002b) 



Irreversible processes are associated with a positive-defined material entropy 



production ( 


Peixoto et al. 


1991 


DeGroot and Mazur 




1984 


Kondepudi and 


Prigogine , 


1998 


Fraedrich and Lunkeit 


2008 


Kleidon 


2009 


). The material 



entropy production is a fundamental quantity in nonequili 



Drium Thermo- 



dynamics (DeGroot and Mazur, 1984 Kondepudi and Prigogine, 1998) and 



gives information about the irreversible processes taking place within the sys- 
tem. General discussions about the entropy budget of the climate system and 
about how to estimate it from climate models can be found in IPeixoto et al.l 



(1991), Goody (2000), Kleidon and Lorenz (2005), Kleidon (2009), Pascale 



et al. (2011a), Pascale et al. (2011b), Lucarini et al. (2011). For a climate 



with a dry atmosphere the material entropy is due to two kinds of processes: 
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dissipation of kinetic energy and sensible heat fluxes, whereas irreversible 
processes associated with the hydrological cycle (Pauluis and Held, 2002af b) 
are not present. If e 2 is the local rate of kinetic energy dissipation such that 
D = f e 2 pdV, the entropy production associated with it reads: 



S, 



kediss 



T 



pdV. 



(14) 



In PlaSim the dissipation of kinetic energy is due to: (i) turbulent stresses in 
the surface boundary layer (which accounts for more than 50% of the overall 
dissipation) and, gravity wave drag, implemented as a Rayleigh friction at 
the highest level with a timescale of 50 days, which we define as D p h ys ; (ii) nu- 
merical dissipation due to numerical diffusion (hyperdiffusion) of momentum 

The total dissipation of kinetic energy 
. Although it is hard to inter- 



(|Johnson[ 1997), which we call D ni 



of the model is therefore D 
pret D num as representative of small scale dissipative processes (jJablonowski 



and Williamson, 2011 ) - the hyperdiffusion schemes do not usually match the 



symmetry requirements of the stress tensor needed to ensure the conservation 



of the angular momentum (Becker, 2001) - these contributions are produced 



by the model and will be taken into account in order to be consistent with 



the model itself (Johnson, 1997 Egger, 1999 Woollings and Thuburn, 2006) 



Sensible heat fluxes are driven by the temperature difference existing 
between the lowermost part of the atmosphere and the surface. The material 
entropy production associated to the sensible heat flux Ft is: 



S„ 



1 



1 



dS, 



(15) 



where T a is the temperature of the first atmospheric level (where Ft is ab- 
sorbed thus heating it) and T$ the surface temperature. The material entropy 
production of the system is therefore: 



The ratio 



'mat 



a 



S, 



kediss 



+ s s 



(16) 
(17) 

which is zero 
dissipation of 

the mechanical energy (Lucarini et al. 2010). The parameter a introduced 



Ssens/ *S kediss 



is a measure of the degree of irreversibility of the system, 
if all the production of entropy is due to the unavoidable 
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above is related to the Bejan number Be as Be = a + 1 (Paoletti et al. , 1989). 
Systems with large a are instead characterized by high thermal dissipation 
relatively to the mechanical viscous dissipation and therefore by a higher 
degree of irreversibility 



4. Results 



4-1. Circulations 

In Fig. Q and @ the cross sections of temperature and zonal wind are 
shown for nine different circulations, namely for TZo = 8 and Tt = 1.5 x 10 -3 
i)}, T s = 10 (Fig. |2(b)fr , F f = 4 x 10 5 (Fig. |2(c)| ) (corresponding 
1/10 and r = 2700s, 3 days, 500 days respectively); TZo = 0.1 and 
(Fig. gdj, J) = 10 3 (Fig. ~ 
and same as before for r); 7?.o = 10 -3 
(Fig. 12(e)), J) = 10 9 (Fig 



(Fig. 
to Q 



10' 



and J 7 / = 



= 10 3 (Fig 
10 (Fig 



2(f) | ) (fi* = 1 
Ijj) = 10 5 
8 and same as before for r) for 



2(f)D (fi* 

the HIGHHC and LOWHC experiments respectively, whereas in Fig. [3] and 
4] the meridional stream function is shown for the same parametric values 
and for HIGHHC and LOWHC respectively. We chose these specific cases 
for illustrative purposes since they correspond to extremal and intermediate 
choices of TZo and Tf. 

Fig. [2(a)| [2(b)] |2(c)l (Fig. [2(j)J gkj) glj] for LOWHC) and Fig. gaj 



2(b) |2 


c) 


(F 


?. 4(a), 


4(b), 



4(c) ) show the low rotation rate case for three 
different intensities of the surface-atmosphere turbulent coupling. Such circu- 
lations are dominated by one Hadley cell in each hemisphere which extends 
northward up to the poles. This is a general consequence of the conser- 
vation of angular momentum and in agreement with the simple theory of 



the Hadley circulation of Held and Hou ( 1980 ) which is able to redistribute 
heat efficiently, 
constant, 



The atmospheric meridional temperature profile is almost 
particularly in the middle atmosphere. This is typical of slowly 



rotating planets (Williams, 1988a Navarra and Boccaletti, 2002), and is due 



to the strong Hadley cell circulation. It is interesting to note the effect of 
the surface drag on shaping the Hadley circulation. For strong surface drag 
(r = 45 min), there is a strong kinetic energy dissipation at the surface, where 
we observe very weak winds (Fig. 2(a), |2(b) ). The intensity of the winds 



increases, and the position of the jet moves further northward as r increases. 
Horizontal shear increases too, and for low surface exchange rate (r = 500 
days) an equatorial westward jet appears. The atmosphere tends to cool 
down since the sensible heat flux from the surface decreases, and therefore 
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the atmosphere is more and more decoupled from the surface - which in turn 
tends to warm up -, and the meridional temperature structure is more and 
more flattened. Comparison of Fig. 3(a) , 3(b) , 3(c) and Fig. 4(a) , 4(b)[4(c) 



reveals that the most vigorous meridional circulation is associated with the 
intermediate case (r = 3 days). Interestingly, for the LOWHC simulations a 
secondary circulation develops within the Hadley cell (Fig. 4(a)). 

In the intermediate range (IZo = 0.1), we have atmospheric circulations 
characterized by strong eastward zonal jets at about 50 — 60° and by a ther- 
mally direct (Hadley) and indirect (Ferrel) meridional cell (Fig. [I] (d,e,f) 
and Fig. [3] (d,e,f)). The general circulation is considerably affected by the 
different surface properties. In particular we note that at large J-f, the 



flow develops strong barotropic horizontal shears, as first discussed by James 



and Gray 



(2(d) 



(1986). Note that none of the three circulations shown in Fig. 
2(f) ) is close to the one we observe on Earth (e.g. Peixoto and Oort 



(1992)). This is not surprising however, since in this simulation setup we 
have removed completely the water vapour, which indeed has a large effect 
on the atmospheric circulation as also shown by these numerical simulations. 
However the climate obtained for IZo pa 0.1, J 7 / ~ 10 3 shows a good simi- 
larity with the mean state of Mars (e.g. see figure 2 of Read (2011)), which 
in fact has a dry atmosphere and relatively similar dimensionless numbers 
(IZo pa 0.2, pa 50). In particular, when we consider the circulation ob- 
tained for Tj pa 50 (Fig. [5]), the resulting patterns are even more similar to 
those of Mars. 

The effect of the surface drag is particularly seen in the meridional cir- 
culation, which is also largely modified by the surface properties. A clear 
thermally direct-indirect cell structure emerges in the intermediate cases 
J-f ~ 10 3 (r ~ 3 days), with the boundaries of the Hadley cell at about 
40°. The intensity and the extent of the indirect cell is greatly reduced in 
the low drag (J 7 / ~ 10 _1 ) case and almost completely suppressed in the high 
drag (J 7 / ~ 10 5 ) case. 

Finally, for the high-rotation runs (TZo = 
tures (Fig. p[g)-(i)), as expected from (10). 
larger contrast in the meridional and vertical profile since the heat transport 
becomes very inefficient, and the thermal structures for high rotation rates 
tend to the radiative equilibrium one. The effect of r is mainly observed in 



10 3 ) we observe small struc- 
The temperature field shows 



the zonal wind profiles 
ional stream function (Fig 




g,h,i) and Fig. [2]^g,h,i)) and in the merid- 
g,h,i) and Fig. [3|g,h,i)). Multi-jet structures 



become more and more evident as the surface drag decreases, as can be seen 
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in Fig. 3(i) and Fig. 4(i 



4-2. Thermodynamic analysis 

In this section we analyze the dissipative properties of the different cir- 



culations described in Sec. |4. 1| as the parameters fl and r, and consequently 
IZo and Tf, are varied. Sensitivity studies of dissipative properties have 



been proposed first by Kunz et al. (2008) and then used extensively in Pas 



cale et al. (2011b) and Boschi et al. (2012) as an insightful way to assess 



the models' tuning and their thermodynamical properties. In the following, 
we plot quantities in the (Q*, r) plane for practical purposes, and we over- 
plot the values of log 10 7?.o and log 10 J 7 /, as can be seen from Fig. 6(a 
Fig. [14(3 



to 



Before starting our analysis, let us note that all steady states are at a 
global mean surface temperature of about 275 K since in a dry atmosphere 
and in a climate without sea ice there are no strong feedback mechanisms 
capable of altering the spatially integrated energy balance at the surface and 
at the top of the atmosphere. 

In Fig. 6(a) and Fig. |7(a)[ ) the total (i.e. due to physical and numer- 
ical processes) dissipation of kinetic energy is shown for the HIGHHC and 
LOWHC experiment. We observe a non-trivial dependence on Q and r. The 
most intense energy cycle is obtained in a region centered around IZo ^ 0.1 
and J 7 / ps 10 3 (r = 3 days and Vl* = 1), with D rj 0.35 WuT 2 . It is interest- 
ing to note how much (~ factor 10) the strength of the general circulation 
is weakened by the lack of the hydrological cycle - the Earth's atmosphere 
has D ~ 2 Wm -2 (Peixoto and Oort, 1992). Simulations of the Snowball 



Earth (Boschi et al. , 2012) have shown that the intensity of the Lorenz en- 
ergy cycle in the "deep frozen" state is about 0.8 Wm~ 2 (differences are 
due to the different vertical resolution and to the fact that the Snowball 
Earth is not completely dry), thus confirming the almost dry nature of the 
Snowball states. Overall, there is a significant decrease of D at low thermal 
Rossby number (TZo < 5 x 10~ 3 ) for any value of JFj and for slow-rotation 
(TZo > 1) and low-drag (J 7 / > 10 4 ) simulations. The dissipation rate remains 
significantly high for slow rotations and for strong drag (J 7 / < 0.1). 

In order to relate the features of D(7Zo, J 7 /) to the characteristics of the 
circulations, we look at the "physical" kinetic energy dissipation D p h ys (see 
discussion in Sect. [3] abo ut the "physical" and "numerical" entropy produc- 
tion) shown in Fig. 6(b) and associated with the boundary layer extraction 

7(b)] ). The two regions 



of kinetic energy from the mean flow (Fig. 6(b) 
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with high surface dissipation are now clearly visible and defined by lines 
of constant J-}, with the first delimited by J 7 / < 0.1 and the second one 
by 10 < Tf < 10 4 . A transition line between the Hadley regime and the 
baroclinic regime is clearly visible for T$ ~ 1. 

The dynamical mechanisms behind these two high-dissipation regions are 
different and related to the very diverse kinds of circulations we have shown in 
Sect. 4.1 For 10~ 2 < TZo < 10" 1 baroclinic waves efficiently convert potential 



energy into kinetic energy (Peixoto and Oort, 1992; Lorenz, 1967 Vallis 



2006) and transport energy polewards. The fact that baroclinic instability 
is strongest for intermediate values of Q is understood from eq. (10): for 
Q — > oo, L c — > - i.e. the size of the baroclinic cyclones becomes so small 
that they are easily dissipated at the surface and do not develop; on the 
other hand, for Q — > 0, L c — > oo and necessarily L c is larger than the 
planet's radius, and therefore the planet cannot accommodate the unstable 
waves. The surface properties have a dramatic impact on the circulation, 
as shown also by James and Gray (1986), because the growth rate of the 
most unstable baroclinic waves is strongly inhibited by horizontal shears 
(James, 1987). Such horizontal barotropic shears build up when there is 
weak surface drag. This explains the drop of D at high Tj and intermediate 
TZo. Conversely, for high surface drag, baroclinicity is also suppressed because 
kinetic energy is rapidly extracted, thus preventing the growth of baroclinic 
eddies. The dissipation (~ 0.16 Wm -2 ) is then associated with the friction 
extracted on the Hadley cell circulation. Our results for 1/5 < Q* < 2, i.e. 
in the baroclinic regime, are in agreement with those of ( |Kleidon et al. 2003 
2006), who indeed considered the case Q* 



1 only. 



In the low rotation regime (TZo > 10) we do not observe local maxima 
in D but monotonic increase with surface drag. This is consistent with the 
fact that in the low rotation regime baroclinicity is mostly absent, and the 
dissipation of kinetic energy is simply related to the strength of the surface 
drag, which extracts kinetic energy from the mean flow, thus causing very 
weak winds near the surface. At high rotation rates (TZo < 10 -3 ) the strength 
of the Lorenz energy cycle tends to zero and becomes insensitive to the surface 
properties. 



Meridional heat transport (Peixoto et al. 1991) is in general a very im 



portant quantity in planetary atmospheres (Lorenz et al. 2001) associated 



with the radiative imbalance between high temperature regions - where the 
emitted infrared radiation is smaller than the absorbed stellar radiation - 
and low temperature regions - where the emitted infrared radiation is larger 
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than the absorbed stellar radiation. The zonal mean of the meridional heat 
transport is worked out by integrating over latitude the longitudinally aver- 
aged TOA radiation budget. A scalar index of the meridional heat transport 
is then defined as half of the sum of the values of the poleward heat transport 
in the two hemispheres at 30° latitude. We also introduce a scalar index for 
the meridional temperature gradient, which in a dry atmosphere is essential 
in controlling the meridional heat transfer, defined as the difference between 
the mean emission temperature Te = (LW/a) 1 ^ of the tropical (30S,30N) 
and the polar (90S, 60S and 60N,90N) regions. The two indices are shown 



in Fig. 8(a) 9(b) 



We observe that meridional heat transport generally increases with IZo, 
in agreement with the results found in Vallis and Farneti (2009). In the 
intermediate-high regime, this result is explained by the inefficiency of the 
baroclinic disturbances in transporting heat, due to their smaller and smaller 
size at large Q* . For Q* > 1/2 the meridional heat transport has a maximum 
for r ~ 3 days. The whole pattern is compatible with the pattern of the 
meridional difference of the emission temperature, which has a minimum in 
the same region (Fig. 8(a) 9(a) ). As the meridional heat transport tends to 
zero at large IZo for any J 7 /, such trends tend to disappear and the meridional 
temperature structure becomes fairly insensitive to the surface properties. 
For slowly rotating simulations (IZo > 10) we still have the largest heat 
transport for r ~ 3 days, although simulations with low r also exhibit large 
heat transports, which may be explained by lower wind velocities in the lower 
branch of the Hadley cell (equatorwards motion). 

In the baroclinic range (1/5 < Q* < 2) the maximum in meridional 
heat transport coincides with the maximum in dissipation/baroclinic activity 
(Fig. 6(a) and 7(a)), whereas in the axisymmetric range (Q* < 1/5) the area 
with the highest dissipation does not correspond with the most vigorous 
meridional heat transport. This is due to the fact that the mechanisms for the 
meridional heat transport are different in the two regimes. In the baroclinic 



range, midlatitude eddies constitute the dominant mechanism (Lorenz ; 1967 



James, 1994). This is particularly clear from the zonal mean of the transient 
eddy flux v'T' (not shown), which reaches the highest values ~ 8 Kms -1 
at 900 hPa and 50 N/S for the values of r maximizing D, compared to 0.5 
Kms" 1 for r = 45 min (at 700 hPa and 60 N/S) and 4 Kms" 1 for r = 500 
days (at 1000 hPa and 50 N/S). Just for the sake of comparison, let us note 

at 850 hPa and 50 N/S 



that for earth's circulation v'T 
(e-g 



x w 15 Kms 

James, 1994). Conversely, in the axisymmetric case, the dominant 
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mechanism of transport is the mean meridional circulation vT through the 
Hadley cell. The meridional stream function (Fig. [3] and [4]) shows that in 
the intermediate case (r = 3 days) the most vigorous meridional circulation 
is observed. 

The Carnot efficiency (Fig. 



10(a) 10(b)) plots reveal that the highest 



values are seen for r ~ 1 — 10 days. We generally observe a certain corre- 
spondence between areas of high efficiency and areas of high meridional heat 
transport, particularly in the intermediate range (1/5 < Q* < 2) where also 
the strongest Lorenz energy cycle is observed. Strong eddies are therefore 
a very efficient way of transporting heat. Likewise, at high rotation rates, 
the drop of meridional heat transport and of the Lorenz energy cycle is also 
accompanied by a strong reduction of the Carnot efficiency. 
Th e entropy production is shown in detail in Fig. 



11(b) 



11 13 Fig. 



show the contribution due to thermal dissipation >S sens (15). This is 



11(a) and 



almost independent of IZo and depends mainly on r. We observe a ridge in 
the values of S sense for r ~ 5 days, which is explained by a trade-off mech- 
anism between the sensible heat flux, which decreases with r independently 
from Q (not shown), and the temperature difference between the surface 
and the near-surface atmosphere, which increases with r also independently 
from Q. The entro py productio n associated with the dissipation of kinetic 
energy, Skediss (Fig. 12(a) 12(b)), depends strongly on both r and Q, and it 



is determined by the pattern of D (Fig. 6 (a" 



7(a)). 



The total material entropy production (16) is the sum of the two, so 
its properties are determined mainly by S sens which is generally larger than 
Skediss (~ 2 — 3 times in the low-medium range, as can be seen in Fig. 14 



where the irreversibility parameter a is shown). The region of maximum 
dissipation («i 5 mWm~ 2 K _1 ) is observed for TZo ~ 0.1 and r ~ 3 — 5 days, 
which corresponds with the steady states with the highest kinetic energy and 
thermal dissipation. Overall, the material entropy production tends to be low 
for high rotation speeds (e.g. IZo ~ 10~ 3 ) and for high r (weak turbulent 
atmosphere-surface coupling) . 

Finally let us note that the difference in the soil heat capacities do not 
significantly alter the patterns of circulation regimes. The most relevant 
difference in the circulation is the splitting of the Hadley cell for low r at 
low Vt* (Fig. [1]). Apart from this we observe an overall weakening of the 
general circulation in HIGHHC case. This is particularly evident in the 
dissipation of kinetic energy (i.e. Lorenz energy cycle, Fig. 6(a) and Fig. 



7(a) ), that in HIGHHC runs is about a factor 0.8 smaller than in LOWHC, 
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in the meridional heat transport (factor 0.9, Fig. 
in the efficiency (factor 0.6, Fig. 10(a) and Fig. 



8(b)] a nd Fig. [9(b)] ) and 
10(b)). The decrease of 



the surface heat capacity is reflected in the decrease of the timescale needed 
to respond to fluctuations of surface energy fluxes, e.g. due to the seasonal 
cycle. In the limit of very small heat capacity the surface temperature would 
adapt instantly to such fluctuations, whereas for infinite heat capacity, the 
surface temperature would not respond to it at all. This loss of positive 
temporal covariance between temperature and diabatic heating is reflected 
in a loss in the generation of available potential energy - which in fact, in 



Lorenz formulation (Lorenz, 1960), is given by the space and time covariance 



between temperature and diabatic heating. 

4.3. Relevance for MEPP 

In this last section we briefly describe our results in the context of the 
MEPP, as this conjecture has gained some momentum in the planetary sci- 
ence community (Lorenz et al. , 2001 Taylor, 2010). In this paper r has been 



considered as a parameter describing the roughness properties of the surface. 
In this respect it is not considered as an internal tunable parameter, as usu- 
ally are other model parameters which are associated with the description of 
subgrid turbulent motions. A different interpretation of r may be found in 
Kleidon et al. (2003) and Kunz et al. (2008), in which r is considered as a 



tunable parameter associated with turbulent processes within the planetary 
boundary layer. For Q* = 1 they showed that the optimum value was the one 
maximizing the material entropy production and the dissipation of kinetic 
energy and claimed this as a strong evidence in support of the MEPP. One 



major criticism against the claim of demonstration of MEPP by Paltridge 



(1975) and Lorenz et al. (2001) was that it did not take into account the 



effects of the rotation speed (Rodgers, 1976; Goody, 2007), which indeed, as 



we have seen, has a large impact on the thermodynamics too. This issue has 



been partially addressed by Jupp and Cox (2010) through a two-box model 



with very simple dynamical assumptions. For sufficiently small, the revised 



model by Jupp and Cox (2010) reproduces to some extent the MEPP states 



obtained without taking the planetary rotation rate into account. 

A well defined maximum in D with respect to r is found in the intermedi- 
ate regime. In this regard, MEPP (when thought as maximum dissipation of 
kinetic energy, as first proposed by [Lorenz (I960)) seems to fit the baroclinic 
regime, i.e. to describe the statistics of the baroclinic macroturbulence. We 
note that the axially symmetric regime, in which the eddy fluxes become 
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unimportant, has no local maximum in the kinetic energy dissipation. Fur- 
thermore, for high rotation rates there is insensitivity to r as far as the 



dissipation of kinetic energy is concerned, in agreement with Jupp and Cox 



(2010) 



A very different answer is obtained when the thermal dissipation is con- 
sidered, which has a well defined maximum with respect to r for any value 
of Q*. Such a maximum arises from the trade-off between the sensible heat 
flux and the surface-atmosphere temperature difference. The total material 
entropy production is mainly determined by the thermal dissipation associ- 



ated with sensible heat flux - which is numerically the largest (Fig. 11, Fig. 



14(a) , Fig. 14(b) ) - and shows its highest values for 1 < r < 10 days. This 



is fairly independent from the rotation speed Q* and determined mainly by 



the surface relaxation timescale r (Fig. 11). In other terms, the possibility 
of applying MEPP and its outcome depend on Q and on the dissipative func- 
tion (e.g. D has no maximum for Q* = 1/10, whereas S mat has a maximum). 
These results confirm the vagueness of MEPP as a physical principle. 

5. Discussion and conclusions 

Stimulated by the ongoing development of exoplanet sciences, in this 
study we have investigated the properties of planetary atmospheres for dif- 
ferent values of two fundamental parameters. The effect of the surface prop- 
erties on the general circulation of a dry Earth-like atmosphere has been 
investigated for different rotation speeds, ranging from Q* = 1/10 to Q* = 8. 



In terms of the dimensionless numbers discussed by Read (2011), this has 
been equivalent to varying the thermal Rossby number TZo from ~ 10~ 3 to 
10 and the friction dimensionless number TZo from ~ 10~ 3 to 10 9 . 

It has been shown that the surface properties, by having a direct effect on 
the turbulent transfer of temperature and momentum between the surface 
and the atmosphere, have a large impact on the dynamics and thermody- 
namics of the general circulation at all the considered rotation speeds. 

For low rotation regimes (TZo > 1) the circulation is essentially axisym- 
metric with no baroclinic eddies. The surface turbulent relaxation timescale 
controls the strength of the meridional circulation, with the most intense 
meridional streamfunction values for J-f ~ 10 (r ~ 3 days). In the medium 
range (10~ 2 < TZo < 1) we have an atmospheric circulation in which baro- 
clinic eddies and the Hadley cell coexist. In this regime r has a dramatic 
effect on the baroclinic eddies, whose formation is inhibited by both very high 
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surface drag - kinetic energy is readily extracted before they can grow - and 
very low surface drag since the building up of strong horizontal barotropic 



shears greatly reduces baroclinicity (James and Gray, 1986 James, 1987) 



This strong effect can also be observed in the meridional circulation, with 
the Ferrel cell disappearing for high and low r. The most vigorous baro- 
clinic activity is observed for 10 < Tf < 10 4 . For fast rotating simulations 
{IZo < 10~ 2 ), Tj controls the multi-jet structure, as we see that the number 
of jet decreases with the increase of the surface drag (Fig. [3] and Fig. [4]). 

The dissipative properties of the different atmospheric circulations show 
very diverse behaviour. The dissipation of kinetic energy, which we diagnose 
directly from the wind field and the associated entropy production, depends 
strongly and non-trivially on both TZo and J 7 / and clearly marks the passage 
from an axisymmetric circulation to a quasi-geostrophic one. The different 
nature of the atmospheric circulations in the axisymmetrix, baroclinic and 
high-rotation rate results in very different trends of the dissipated kinetic 
energy. In particular we see that in the axisymmetric slow-rotating case the 
dissipation of kinetic energy increases monotonically with the surface drag. 
Conversely, in the intermediate baroclinic range the energy cycle peaks for 
10 2 < J 7 / < 10 3 where a clear peak is observed. Such trends persist at very 
low IZo (< 10~ 2 ), but the strength of the atmospheric energetics is greatly 
reduced so that the peak tends to flatten completely for IZo > 10 2 . 

The thermal dissipation S sens is instead fairly insensitive to IZo and is 
determined mainly by the timeconstant r, due to a trade-off mechanism 
between the temperature difference and the heat flux. A very interesting 
pattern is found for the meridional heat transport and the Carnot efficiency. 
In the intermediate-high rotation regimes, circulations with the largest dis- 
sipation of kinetic energy have also the highest efficiency and the largest 
meridional heat transport. This result highlights the importance of the baro- 
clinic eddies in the determination of these main thermodynamic quantities. 
In the low rotation regimes instead, high meridional heat transport does not 
necessarily imply high dissipation of kinetic energy or the highest efficiency. 

Let us also note the effect of the two different kinds of surfaces we have 
considered: HIGHHC, with ocean and albedo heat capacity, and LOWHC, 
with albedo and heat capacity typical of rock. In spite of similar dissipative 
properties, we observe that LOWHC has generally a more intense circulation 
and larger Carnot efficiency. The weakening of the circulation is associated 
with the loss of positive temporal covariance between temperature and di- 
abatic heating, which in turn is reflected into a loss in the generation of 
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available potential energy. Also a curious feature in the circulation is the 
splitting of the Hadley cell for low r at low Q* (Fig. [4]). 
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Q* TZo max(J-j) min(J-/) 

1/10 7.98 1.5 x 10~ 3 3.9 x 10 5 

1/5 1.99 6.2 x lO" 3 1.5 x 10 6 

1/2 0.31 3.8 x lO" 2 9.9 x 10 6 

1 0.08 0.15 3.9 x 10 7 

2 1.9 x lO- 2 0.62 1.6 x 10 8 
4 4.9 x 10~ 3 2.5 6.3 x 10 8 
8 1.2 x 10~ 3 9.9 2.5 x 10 9 



Table 1: Range of variation for TZo and Tj. For each value of the rotation rate the 
minimum and maximum value of J 7 / (corresponding to the smallest and largest value of r, 
see definitions Q and ^) are shown. A value of A6h — &9hE ~ 60 K has been chosen. 
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ures' captions 
» Figure [T] 

Zonal mean of temperature and zonal wind (HIGHHC) for TZo = 8 
(fl* = 1/10) (r = 2700s (a), 3 days (b), 500 days (c)), fT = 1 (r = 
2700s (d), 3 days (e), 500 days (f)), O* = 8 (r = 2700s (g), 3 days (h), 
500 days (i)) 

• Figure [2] 

Zonal mean of temperature and zonal wind (LOWHC) for fl* = 1/10 
(r = 2700s (a), 3 days (b), 500 days (c)), fi* = 1 (r = 2700s (d), 3 
days (e), 500 days (f)), Q* = 8 (r = 2700s (g), 3 days (h), 500 days 

(i)) 
» Figure [3] 

Meridional mass transport streamfunction (10 9 Kgs^ 1 ) of the HIGHHC 
experiment for Q* = 1/10 (r = 2700s (a), 3 days (b), 500 days (c)), 
Q* = 1 ( r = 2700s (d), 3 days (e), 500 days (f)), fi* = 8 (r = 2700s 
(g), 3 days (h), 500 days (i)). 

• Figure [4] 

Meridional mass transport streamfunction (10 9 Kgs -1 ) of the LOWHC 
experiment for Q* = 1/10 (r = 2700s (a), 3 days (b), 500 days (c)), 
n* = 1 (r = 2700s (d), 3 days (e), 500 days (f)), tt* = 8 (r = 2700s 
(g), 3 days (h), 500 days (i)). 

• Figure [5] 

(a) Zonal mean of temperature and zonal wind and (b) meridional mass 
transport stream function (10 9 Kgs -1 ) for TZo ~ 0.1, J 7 / ~ 10 3 show a 



good similarity with the mean state of Mars (e.g. see figure 2 of Read 
( 201 ![ )), which in fact has a dry atmosphere and close dimensionless 
numbers (TZo w 0.2, Tf w 50). 

Figure |6] 

Total kinetic energy dissipation for (a) HIGHHC and (b) LOWHC 
experiments; contribution to the total kinetic energy dissipation due 
to parametrizations representing physical processes (boundary layer 
stresses and gravity wave drag) for the HIGHHC (c) and LOWHC (d) 
experiment. Overp lotted (as in all the following plots) are the values 
of log 10 7?.o (dotted) and log 10 J-/ (continuous). 
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Figure [8] 

Equator-to-pole temperature difference (a) and atmospheric meridional 
energy transport (b) for the HIGHHC simulation. 

Figure [9] 

Equator-to-pole temperature difference (a) and atmospheric meridional 
energy transport (b) for the LOWHC simulations. 



Figure 10 



Carnot-like efficiency for the HIGHHC (a) and LOWHC (b) experi- 
ments. 



Figure [IT] 

Entropy production associated with surface sensible heat flux (HIGHHC 
(a) and LOWHC (b)). Units in 10" 3 Wm^K" 1 , 



Figure [12] 

Dissipation of kinetic energy (HIGHHC (a) and LOWHC (b)). Units 
in 10- 3 Wm-'r 1 . 



Figure 13 



Total material entropy production (HIGHHC (a) and LOWHC (b)). 
Units in 10" 3 Wm^K" 1 , 



Figure [14] 

Irreversibility parameter a for (a) HIGHHC and LOWHC experiments. 
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Figure 1: 






(a) T f = 1.5 x 10" 3 , TZo = 8 (b) F f = 10, TZo = 




(c) 



(d) Ff = 10~ l , TZo = 10 _1 (e) J 7 / = 10, TZo = 10 (f) J/ = 4 x 10 5 , TZo = 10 






J/=4x 10 5 , TZo 
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(g) J 7 / = 10, TZo = 10" 3 (h) F f = 10 5 , "Ko = 10~ 3 (i) F f = 10 9 , TZo = 10~ 3 






(j) Tf = 1.5 x 1CT 3 , IZo = 8 (k) T f = 10, Tlo = 8 (1) J/ = 4 x 10 5 , Ko = 
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(m) F f = 1CT 1 , Ko= 10" 1 (n) J 7 / = 10 , 7£o = 10 (o) F f = 4 x 10 5 , Tlo = 10" 






(p) J 7 / = 10, fto = 10~ 3 



(q) J) = 10| 4 Tlo = 10- 3 
Figure 2: 



(r) F f = 10 9 , fto = 10~ 3 
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(g) J> = 10, TZo = 10~ 3 



(h) Tf = lOj^ TZo = 10~ 3 (i) Tf = 10 9 , TZo = 10~ 3 
Figure 3: 
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(g) J> = 10, Ko = 10- 3 



(h) Tf = lOj^g Ko = 10- 3 (i) Tf = 10 9 , Ko = 10~ 3 
Figure 4: 
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KE dissipation [ 10" 1 W m" 2 ] - LOWHC 
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Emission temp, eq-pole diff [K] - HIGHHC 
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Emission temp, eq-pole diff. [K] - LOWHC 
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Carnot efficiency [%] - HIGHHC 
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sensible heat flux EP [10 3 W m" 2 K" 1 ] - HIGHHC 
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KE dissipation EP [1 0" 3 W m" 2 K" 1 ] - HIGHHC 
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IRREVERSIBILITY PARAMETER - HIGHHC 
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